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Multiple-frequency and multiple station radio wave absorption data from the auroral zone are 
used to investigate the electron density content in the D region during a certain type of absorption 
events. The type is characterized by smooth variability and widespread geographic extent; the ap- 
proximate electron density profiles are derived and discussed. The general features of the profiles 
were found compatible with the ionization due to the high energy electrons and the associated 
bremsstrahlung x-rays. Such events ar^^ also discussed in relation to the visual auroral and mag- 
netic records. 



1. Introduction 

One of the more significant questions engaging the 
attention of the radio-physicists in the auroral zone 
concerns the depth and extent of the radio wave ab- 
sorbing region. Among the recent outgrowths of the 
earlier single frequency observations (about 30 Mc/s) 
by Little and Leinbach [1958], and Chapman and Little 
[1957] is the introduction of the multiple-frequency 
technique jointly by the CRPL and the Geophysical 
Institute [Parthasarathy, Lerfald, and Little 1963; 
Little, Lerfald, and Parthasarathy, 1964]. This paper 
will seek to present some of the results of the multiple- 
frequency and muhiple-station program, relevant to 
the associated radio wave absorption during the 
dawn-breakup phase of the auroral displays. 

2. Method of Deriving Ne{h) Profiles 

The method of deriving the ionization profile was 
discussed in some detail by Parthasarathy et al. [1963], 
and will only be outlined and supplemented here. The 
radio w^ave absorption from any ionized stratum is 
a function of the wave frequency, the electron gyro 
frequency, and in the case of the nondeviative absorp- 
tion, is also proportional to the electron density. For 
a very narrow beam, zenith directed antenna. 
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where the absorption as measured by the antenna is 
A^ the K function is absorption per unit path length 
when the electron density is unity, and A^ is the electron 
density at height h. 



The specific absorption K is a. function of the elec- 
tron collision frequency, i^m (appropriate to the electrons 
with the most probable velocity), the wave frequency, 
/, the angle of propagation, (/), with respect to the mag- 
netic field, and the electron gyrofrequency, s. The 
form of this function used in the present work is due 
to Sen and Wyller [1960], which takes into account the 
dependence of the collision frequency of each electron 
on its energy [Pack and Phelps, 1961; Phelps I960]. 
The function shown in figure 1, for a set of radio fre- 
quencies, is from a tabulation by Jean Chorbajian 
et al. [1962], by considering the tabulation merely as 
one of K(f, Vm, 5, 0). The collision frequency Vm[h) 
required to convert K{f, Vm, 5, 0) into K{f, h) was ob- 
tained by combining the momentum transfer cross 
section of Pack and Phelps and the atmospheric gas 
pressure obtained by rocket data [see, Nordburg, 
1961]. The simple relationship is i/,„ =6.3 X lO^P, 
where P is gas pressure in dynes per cm^. 

Under certain favorable conditions revealed by the 
frequency dependence of absorption, a given N(h) 
function may be tested through the integral relation- 
ship (1) to see whether it satisfied the observed data 
in those frequencies. The larger the frequency range 
of simultaneous data that is available, the more con- 
fidence can one place in the detailed validity of the 
profile. As regards the type of N{h) profiles that need 
be considered, the following observational and theo- 
retical surmises may be pointed out. The differential 
energy spectrum of the downflowing flux, whether 
protons or electrons, is expressible by a simple power, 
or exponential law; the atmospheric gas density as 
well as the effective loss coefficient of electrons by 
eventual recombination monotonically decrease with 
height. Under the circumstances, the resulting elec- 
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tron density profile may be reasonably expected to 
approximate a low degree polynomial of height. (That 
a low degree polynomial is a satisfactory representa- 
tion of the actual ionization profile may be seen by 
reference to the proton-derived profiles of Webber and 
Freier [1963] and the electron- and x-ray-derived pro- 
files of Aikin and Maier [1963].) Basically, the pro- 
file derivation outlined in Parthasarathy et al [1963], is 
only an application of this; the profile was derived by 
assuming it to be a low degree polynomial (cubic) and 
then computing the coefficients of the polynomial that 
could best fit the absorption data. 

A feature of the auroral latitudes is the seasonal 
variation of the gas density and pressure in the D 
region, which is shown in figure 1. Because of the 
trivial variabihty of s and </) with height in the zenith 
D region, it is easy to revise the K{h) obtained from a 
given Pm{h), in order to account for the seasonal 
variation in Pm{h). 



3. Experimental Details 

Cosmic noise level was monitored at College, Alaska, 
in eight effective frequencies of 5~, 5+, 10~, 10^, 20~, 
20^, 30, and 50 Mc/s, where the + and — signs desig- 
nate the ordinary and extraordinary modes, respec- 
tively. The linearly polarized antennas consisted of 
a pair of 0.4X dipoles spaced 0.4X apart and 0.15X 
above a reflecting mesh of 1.25k^. The beam thus had 
rotational symmetry about the zenith of 60° between 
half power points. The circularly polarized antennas 
had the form of a square consisting of four 0.4X dipoles, 
the other details being the same as the linear system. 
Conventional ring circuits enabled the channeling of 
the right circular and left circular noise intensities into 
separate channels, with the appropriate minimum level 
detectors to reject the atmospherics. For additional 
experimental details, reference may be made to 
Parthasarathy and Lerfald [1965]. 
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Figure 1. Upper atmospheric pressure and density for the summer and winter seasons (S and W) for the auroral zone, after CIRA, 1961. 

The electron-neutral gas collision frequency p,,, for the two seasons was obtained by combining the gas pressure data and the momentum transfer cross section of Hack and Phelps 
[1961]. The specific absorption Kif, h) is appropriate to the summer conditions; shifting the curves downwards by 5 km gives a close enough approximation for the winter K function. 
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The multiple station network consisted of 27.6 Mc/s 
cosmic noise receivers at College, Ft. Yukon, Kotze- 
bue, and Healy, Alaska. 

4. Data Analysis 

The quiet day levels for the dual polarized frequen- 
cies were obtained as detailed by Little et al., [1964]. 
Essentially, it consists of deriving the absorption-free 
noise level at any sidereal hour from the fact that, if 
the absorption were zero, the two circular modes would 
register equal levels. This dual polarization method of 
deriving the true level may be considered as one of the 
more promising techniques available for HF absorp- 
tion studies. The conventional technique of plotting 
the levels (from a full year of data) against sidereal 
time and then drawing the top envelope was used in 
the case of the linearly polarized frequencies. The 
raw, wide beam absorption data were then reduced to 
the equivalent pencil beam data. The reduction is 
straightforward, providing the distribution of ioniza- 
tion over the major part of the beam could be verified 
to be approximately uniform. 

Auroral absorption events obtained from the Alaskan 
network of stations were recently analyzed by Berkey 
and Parthasarathy [1964] who found it possible to 
group them into a few distinct types. Events which 
appear isolated and show smooth and slow variabihty 
were designated as Type S events. Typical examples 
may be seen in figures 3a and 3c. They occur mostly 



during predawn hours (fig. 2) and invariably in the 
midst of a few hours of little or no absorption. No 
seasonal trend was evident during at least the period 
1961-1963. They were decidedly more common in 
the southern boundary (College) of the auroral zone 
than in the midauroral station of Ft. Yukon. 

Comparison with the magnetograms showed that 
these Type S events could not always be associated 
with a magnetic bay; and whenever they were asso- 
ciated, it was found that the magnetic bay always 
preceded by several minutes. There was poor cor- 
relation between the magnitudes of the absorption 
event and the magnetic bay. 

Comparison with the all-sky camera pictures showed 
that the Type S events occurred mostly during the 
breakup phase of the arcs in the dawn meridians. 
Typically, the arc becomes progressively brighter, 
moves slightly equatorward and then breaks up without 
violent motion; the resulting cloudlike patches, to- 
gether with the striking feature of a dark, well-bounded 
fl-shaped band near the northern side, then move 
eastward. This phase of the auroral display has been 
investigated in more detail by Akasofu [1964]. Of 
relevance to the mechanism of precipitation would be 
the fact that the histogram of the number of cases 
versus magnitude of absorption (in steps of 0.5 dB) 
showed a distinct peak at around 1.5 to 2.0 dB at 30 
Mc/s. For greater details of the analytical results 
reference may be made to Berkey and Parthasarathy 
[1964]. 
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Figure 2. Histogram of occurrence frequency of Type S events in 
the four Alaskan stations during 1961-1963. 

The observations al (lollej^e and F't. Yukon extended to 30 months; Kotzebue, 18 months; 
and Healy, 9 months. 
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Figure 3. The multiple frequency charts corresponding to the samples of derived electron density profiles shown in Figure 6. 

Note the typical shielding action of the F layer during the daytime leading to serious apparent absorption in the lowest of the frequencies, 5 — and 5 + Mc/s. Absorption above about 
30 dB cannot be accurately measured in these frequencies. Within the dynamic range, the charts are dB-linear, except the 30 and 50 Mc/s channels, which are power linear. Obvious 
symptoms of interference may also be recognized in some of the charts. 



Figure 3 shows several examples of these auroral 
absorption events as registered in the eight effective 
frequencies. The geographical extent of the events 
and hence the approximate uniformity of ionization 
may be seen in figure 4 which shows the 27.6 Mc/s 
absorption in College, Ft. Yukon, and Kotzebue, the 
latter two being approximately 200 km (north) and 700 
km (west) respectively. No data from stations south 
of College were available for the events in figure 4, 
but whenever data for similar events were available 
from Healy (100 km to the south), only minor differ- 
ences would be noticed between the absorption at 
College and at Healy. 

The reduced pencil beam data for a few instants 
are plotted in figure 5. Using the values of absorp- 
tion read on the best fitting line for the six or eight 
effective frequencies (rather than the acutal reduced 
data), the instantaneous ionization profiles shown in 
figure 6 were then derived. It may be emphasized 



[Paratharathy et al., 1963] that the profiles tend to be 
somewhat uncertain beyond about 80 km. However, 
by multiplying the N{h) profiles with the winter K- 
functions for any of the frequencies between 5 and 50 
Mc/s, it may be seen that the observed absorption datum 
in each of the frequencies is fully (> 95%) account- 
able by the profiles in the region below the 85 km 
stratum. Arbitrary models of alternate profiles which 
do not ensure a dominant absorbing region below 85 
km would be incapable of yielding the observed fre- 
quency dependence of absorption in figure 5. The 
exact shape of the profile in the region contributing 
only a minor fraction to the data (>80 km) is open 
to arbitrary adjustments, subject only to the require- 
ments of physically plausible shapes, and the inte- 
grated absorption in the adjustable region remains 
small. Because of their sensitivity to regions 
above 80 km, the need for data below 5 Mc/s to sup- 
plement the HF data can hardly be overemphasized 
in ensuring more accurate higher altitude profiles. 
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Figure 4. The 27.6 Mcjs records from the other Alaskan stations of Ft. Yukon (F), College (C), Kotzebue (30 Mcjs) (K), and the conjugate 
area station in Macquarie Island (m); (a), (b), and (c) are during the times of the corresponding ones in figure 3. 



5. Discussion of Results 

A few obvious inferences may be made from the 
profiles, appropriate to these isolated, smooth, and 
geographically widespread absorption events. Pro- 
files which extend down to these low altitudes seem 
compatible with the ionization due to the electrons 
and the associated x-ray Bremsstrahlung, such as 
those proposed by Chapman and Little [1957] and 
Aikin and Maier [1963]. Ansari [1963] has pointed 
out that electrons of 10 keV energy are inefficient 
sources of x-rays; detailed calculation of x-ray pro- 
duction rates has also been made by Rees [1964]. 
It was not found possible, however, to propose from 
among these (and the many other computed profiles) 
a 'model' profile. A statistical analysis, using a lim- 
ited amount of data (from 25 Type S events), showed 
that the events may be defined in the frequency range 
of 5 to 50 Mc/s by the relation k(f)=K'f-^. The 
values of n ranged between 1.3 and 1.9, with a mean of 
1.5. This mean value may be contrasted with the 
value of 1.8 when absorption events were considered, 
without classifying them into the several types [Lerfald 
et al., 1964] . The obvious inference can thus be made 



that during the dawn-breakup phase of the aurora, 
the primary particle spectrum is typically much harder 
than at other times. Of major significance would be 
that during any auroral display, and decidedly so during 
the breakup phase, the region relevant to HF-VHF 
absorption is well below the visual auroral altitude of 
about 100 km and extends down to the 50 km region. 
An essentially similar inference was earlier made by 
Little and Leinbach [1958] by comparing the /-min 
data and the simultaneous cosmic noise absorption 
data at 30 Mc/s. 

Parthasarathy and Berkey [1965] had also found that 
the frequency exponent during the Type F events was, 
on an average, again about 1.5. The Type F events 
consist of those with a fast, flash-like onset of absorp- 
tion and were found to be typical of the midnight 
breakup phase of quiet auroral forms. 

The Ne{h) profiles may be used to derive the energy 
spectrum of the primary electrons if the rate coeffic- 
ients of recombination, attachment, etc., are known as 
a function of height. The effective recombination 
coefficient, a^^^h), may be derived by using the basic 
laboratory data in conjunction with the atmospheric 
parameters [see, for example, Webber, 1962] . The 
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Figure 5. The reduced pencil beam absorption versus frequency. 
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Figure 6. Electron density profiles derived by using the observed 
radiowave absorption data from figure 5 and the winter param- 
eters in figure 1 (see text). 

alternative approaches (to be published) would be to 
derive the coefficient by two independent methods: 
(a) By means of a set of polar cap event (PCE) electron 
density profiles, and the corresponding sateUite meas- 
ured proton flux data; and (b) by means of the Ne(h) pro- 
files obtained during the recovery phase of a special 
category of absorption events characterized by their 
flash-like onset and rapid self-recovery. 



The multiple frequency program and its apphcation 
to profile derivation are due to the very close collabo- 
ration between one of us (R. P) and C. G. Little and 
G. M. Lerfald of the CRPL. We acknowledge the 
enthusiastic help of Lillian Bogert and D. Ruff in the 
scaling and data reduction processes. The investi- 
gation was financed partly by the Advanced Research 
Projects Agency under Contract No. 186-61 and partly 
by the National Science Foundation under Grant 
GP-947. 
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